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(_^ 1 Abstract. In this paper we enumerate fc-noncrossing RNA pseudoknot structures with given 

— ^ ' minimum arc- and stack-length. That is, we study the numbers of RNA pseudoknot structures 

with arc-length > 3, stack-length > u and in which there are at most fc — 1 mutually crossing 

[31 
bonds, denoted by Tj.^(n). In particular wc prove that the numbers of 3, 4 and 5-noncrossing 

RNA structures with arc-length > 3 and stack-length > 2 satisfy TgjC'T') ~ iCs n~^2.5723", 

PQ ! T^^U") ~ ^^"4™"^ 3.0306", and T!;^l,(n) ~ X5 n-^S 3.4092", respectively where K^^.K^.K^ 

Q ' are constants. Our results arc of importance for prediction algorithms for RNA pseudoknot 

* ^ , structures. 

(N 

K^ ' 1. Introduction 

CN ' An RNA structure is the helical configuration of an RNA sequence, described by its primary se- 

quence of nucleotides A, G, U and C together with the Watson-Crick (A-U, G-C) and (U-G) 

1^^^ ■ base pairing rules specifying which pairs of nucleotides can potentially form bonds. Subject to 

C^ . these single stranded RNA form helical structures. The function of many RNA sequences is often- 

times tantamount to their structures. Therefore it is of central importance to understand RNA 

y^ ' structure in the context of studying the function of biological RNA, and in the design process of 

}_( ' artificial RNA. In this paper we enumerate fc-noncrossing RNA structures with arc-length > 3 and 

stack-length > ct, where ct > 2. The main idea is to consider a certain subset of fc-noncrossing 
core-structures, that is structures with minimum arc length 2, in which there exists no two arcs of 
the form (i, j), (i-|-l, j — 1) and no arcs of the form (i,i-|-2) with isolated i-l-1. We prove a bijcction 
between this subset of core structures with multiplicities and fc-noncrossing RNA structures with 
arc-length > 3 and stack-length > cr, where cr > 2. Subsequently, we derive several functional 
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equations of generating functions, based on which transfer theorems imply our asymptotic formu- 
las. The paper is relevant for prediction algorithms of pseudoknot RNA since it proves that the 
numbers of fc-noncrossing RNA structures with arc-length > 3 and stack-length > a exhibit small 
exponential growth rates. The results suggest a novel strategy for RNA pseudoknot prediction. 



2. Diagrams, matchings and structres 

A diagram is labeled graph over the vertex set [n] = {1, ... ,71} with degree < 1, represented by 
drawing its vertices l,...,n in a horizontal line and its arcs (i,j), where i < j, in the upper 
halfplane. The vertices and arcs correspond to nucleotides and Watson-Crick (A-U, G-C) and 
(U-G) base pairs, respectively. 





10 11 12 13 14 15 



Figure 1. fc-noncrossing diagrams. Top: 3-noncrossing diagram with arc-length > 3, 
(2,5), (7,10), the arc (7,10) being isolated. Hence we have a 3-noncrossing diagram 
A = 3, cr = 1 diagram without isolated vertices. Bottom: 3-noncrossing (no red/purple 
cross), X — 4, a — 2 diagram with isolated vertices 13. 



We categorize diagrams according to the maximum number of mutually crossing arcs, fc — 1, the 
minimum arc-length. A, and the minimum stack-length, a. Here the length of an arc (i,j) is j — i 
and a stack of length <t is a sequence of "parallel" arcs of the form ((i,j), (i + 1, j — 1), . . . , (i + (cr — 
1), j — (cr — 1))). In the following we call a fc-noncrossing diagram with arc- length > 2 and stack- 
length > cr a fc-noncrossing RNA structure. We denote the set (number) of fc-noncrossing RNA 



RNA PSEUDOKNOT STRUCTURES WITH ARC-LENGTH > 3 AND STACK-LENGHT > o- 



structures with stack-size > cr by Tk_a{n) (Tj,g.(n)) and refer to fc-noncrossing RNA structures 
for fc > 3 as pseudoknot RNA structures. A /c-noncrossing core-structure is a fc-noncrossing RNA 
structures in which there exists no two arcs of the form {i,j),{i + l,j — l). We denote the set 
(number) of core-structures having h arcs by Ck{n,h) {Ck{n,h)) and Ck{n) {Ck{n)) denotes the 
set (number) of core-structures. The set (number) of RNA structures with arc-length > 3, is 

[31 [31 

denoted by T^ ^(n) (TJ. ^(n)). For fc = 2 we have RNA structures with no 2 crossing arcs, i.e. the 
well-known RNA secondary structures, whose combinatorics was pioneered by Waterman et.al. 
[IllliniEOlESlEl]. RNA secondary structures are r2,i(n)-structures. We denote by fk{n,l) the 
number of fc-noncrossing diagrams with arbitrary arc-length and (. isolated points over n vertices. 
In Figure [2] we display the various types of diagrams involved. 






Figure 2. Basic diagram types: (a) perfect matching (/3(8,0)), (b) partial matching 
with 1-arc (5, 6) and isolated points 2, 7 (faiS, 2)), (c) structure with arc-length > 3 and 

[31 

stack-length > 2 and no isolated points (Tg 2(8)) and (d) structure with arc- length > 3, 
stack-length > 3 and isolated points 4, 5 (Tj 3(8)). 



The following identities are due to Grabincr et.al. [5] 



(2.1) 



(2.2) 




det[V,(2x)-/,+,(2x)]t^7ii 



e"" det[I,-j{2x) - I^+j{2x) 



ifc-i 



where Ir{2x) = J2j>o jUr+i)] denotes the hyperbolic Bessel function of the first kind of order 
r. Eq. (|2.ip and (|2.2p allow "in principle" for explicit computation of the numbers fk{n,€). In 
particular for fc = 2 and fc = 3 we have the formulas 



(2.3) /2(n,^)= ( jq„-«)/2 and hi^A) ^ 



Cr^^^Cr, 



CU 



+1 
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where Cm denotes the 77i-th Catalan number. f^{n,t} resuhs from a determinant formula enu- 
merating pairs of nonintersecting Dyck-paths. In view of fk{n, i) = {^i) fk{n — i, 0) everything can 
be reduced to perfect matchings, where we have the following situation: there exists an asymp- 
totic approximation of the hyperbolic Bessel function due to |17j and employing the subtraction 
of singularities-principle [T7j one can prove 

(2.4) VfceN; h{2n,0)^^k{n)(- 

where pk is the dominant real singularity of X]n>o fk{'2n,0)z"' and ipk{n) is a polynomial over n. 
Via Hadamard's formula, pk can be expressed as 

(2.5) pk= lim(/fe(2n,0))-^ . 

n — >QO 

Eq. (|2.4p allows for any given fc > 2 to obtain ipk (n) , explicitly. 

As for the generating function and asymptotics of A:-noncrossing RNA structures we have the 
following results from [lOl [TTJ [12] . First the number of fc-noncrossing RNA structures with (^^) 
arcs, Tfc^i(n, ^^), and the number of fc-noncrossing RNA structures, Tk,i{n), are given by 

(2.6) Tfc,i(n,^) = ^(-l)''("^^)A.(n-26,£) 

L"/2J / _ /,\ f "~26 ^ 

(2.7) T,,i(n) = E(-1)V5 E^(^-26,£) , 

where {J2l=o Ik{n — 25, £)} is given via eq. (|2.2p . Secondly we have 

10.4724.4! /^ 5 + V2T V 

('■') ^^^^^") ^ n(n-l)...(n-4) \~^^l 

(2.9) J^\in) ^ ^lilllLiL^ 4.54920" . 

n[n — 1) ... (71 — 4) 



The particular class of fc-noncrossing core-structures, i.e. structures in which there exists no two 
arcs of the form (i, j), (i -|- 1, j — 1) will play a central role in the following enumerations: 

Theorem 1. [T3] (Core-structures) Suppose k G N, k > 2, let x be an indeterminant, p^ the 
dominant, positive real singularity o/ ^^^^^ /fc(2ri, 0)z" (eq. i2.5\) ) and ui{x) ~ . ^ .2 ■ Then the 
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(2.10) Q(n,/i) = ^(-l)''-''-if''/)T,^i(n-2;i + 26 + 2,5+l 



numbers of k-noncrossing core- structures over [n], Ck{n) are given by 

/i-i 

fc=0 ^ 

Furthermore we have the functional equation 

(2.11) EC.N-" - , ^ V/,(2n,0)f f'^T 

^ — ' U\X^ — X + 1 ^ — ' \U\X^ — X + 1 I 

n>0 ^ n>0 ^ ' 

and 

(2.12) C,(n)~^,(n) (— 

where Kk is a dominant singularity o/ ^^^^^ Cj,(n) and the minimal real solution of the equa- 
^^'^^ u x"-x+i ~ P'^ '^^^ '^^ ("") *'^ '^ polynomial over n derived from the asymptotic expression of 
fk[2n,Q) ^ ^n) (£)'' of eq. (2^. 



The following functional identity [TT] relates the bivariate generating function for Tk.i{n,h), the 
number of RNA pseudoknot structures with h arcs to the generating function of fc-noncrossing 
perfect matchings. It will be instrumental for the proof of Theorem |4] in Section H) 

Lemma 1. Let fc e N, fc > 2 and z, u be indeterminants over C. Then we have 

n>Qh<n/2 n>0 ^ ' 

In particular we have for u = 1, 

1 / \2n 

(2^4) j;t,,(„),-.-,— ^5:a(2„,o) -,-i^ , 

n>0 n>0 ^ ' 

In view of Lemma [1] it is of interest to deduce relations between the coefficients from the equality 
of generating functions. The class of theorems that deal with this deduction are called transfer- 
theorems [7]. One key ingredient in this framework is a specific domain in which the functions 
in question are analytic, which is "slightly" bigger than their respective radius of convergence. It 
is tailored for extracting the coefficients via Cauchy's integral formula: given two numbers 0, R, 
where R> I and Q < cj) < ^ and p e M the open domain Ap{(f), R) is defined as 

(2.15) Ap(0, R)^{z\ \z\ <R,z^p, |Arg(z - p)\ > 0} 
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A domain is a Ap-domain if it is of the form Ap((/), R) for some R and (f>. A function is Ap-analytic 
if it is analytic in some Ap-domain. We use the notation 

(2.16) {f{z) ~ O {g{z)) as z ^ p) <=^ [f{z)/g[z) is bounded as z — > p) 

and if we write f{z) — 0{g{z)) it is implicitly assumed that z tends to a (unique) singularity. 
[z"] f{z) denotes the coefficient of z" in the power series expansion of f{z) around 0. 



./W = E 



n>0 ""^ 



Theorem 2. [H] Let f{z), g{z) he a Ap-analytic functions given by power st 

and g{z) = X]n>o^"-^"- Suppose f{z) = 0{g{z)) for all z G Ap and bn ^ ip{n){p^^)"' , where (p{n) 

is a polynomial over n. Then 



(2.17) 

for some constant K . 



[z"]/(z) ~ K [z-]g{z) = Kbn^ K^{n){p-'Y 



3. Exact Enumeration 

Our ffi'st result, Theorem [3l enumerates /s-noncrossing RNA structures with arc-length > 3 and 
stack-length > a. The structure of the formula is the exact analogue of the Mobius inversion of 
eq. (|2.10p [T5] . which relates the numbers of all structures and the numbers of core-structures: 



Tfc.aC'^'M =EbJ-i 



(b+{2-a)(h-b)- 
\ h-b-1 



)Cfc(ri — 2&, h— b). While the latter cannot be used in order 



to enumerate A;-noncrossing structures with arc-length > 3, see Figure [S] the set 




i-1 i i+3 i+4 

\*— distance=4- 




|-*distance=2*] 



Figure 3. Core-structures will in general have 2-arcs: the structure 5 £ r3_2(12) (Ihs) 
is mapped into its core c{S) (rhs). Clearly S has arc-length > 4 and as a consequence of 
the collapse of the stack {{i + l,j + 3),{i + 2, j -|- 2), (i -I- 3, j)) (the purple arcs are being 
removed) into the arc {i + 3, j), c{S) contains the arc (i, i + 4), which is, after relabeling, 
a 2- arc. 
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(3.1) C^(?i, /i) = {(5 I (5 G Cfe(n, /i); /] (i,i + 2); i + 1 is an isolated vertex. } 
turns out to be the key. 

Theorem 3. Suppose we have k,h,(7 G N, /s > 2, h < n/2 and a > 2. Then the numbers of 
k-noncrossing RNA structures with arc-length > 3 and stack-length > a having h arcs is given by 

(3.2) Tgl(.,/.)= E (^+(\-:f-^)-^)Q(n-2M-6) 
where C'^{n,h) satisfies 

(3.3) Q(n, h) = Y,{-iy ('' ~ ^^) Cfe(n - 3j, h - j) 

j= 

and Ck{n',h') is given by Theoreni\^ 



^ ' 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


Ti'i(n) 


1 


1 


1 


1 


1 


2 


4 


9 


19 


40 


82 


167 


334 


682 


1398 


2917 


6142 


13025 


Ti^i('^) 


1 


1 


1 


1 


1 


1 


1 


2 


4 


8 


14 


24 


40 


68 


118 


209 


371 


654 



Proof. We consider C^ (n, h) (eq. (|3.ip ) and call an arc (i, i + 2) with isolated i + 1 a bad arc. It is 
straightforward to show that there are ("~ •') ways to select j bad arcs over [n\. Since removing a 
bad arc by construction removes 3 vertices we observe that the number of configurations of at least 
j bad arcs is given by ("~ -') C/c (n — 3j, h — j). Via the inclusion-exclusion principle we accordingly 
arrive at 

(3.4) q(n,M = E(-l)'(""-^0cfe(n-3j,/i-j) . 

We next observe that there exists a mapping from /c-noncrossing structures with h arcs with arc- 
length > 3 and stack-length a >2 over [n] into UCT-i<b</i-i^fe('^ — 2b,h — b): 

(3.5) c:Tf}^{n,h)^[}^^^^^_Cl{n~2b,h~b), 5 ^ c{5) 

which is obtained in two steps: first induce c{5) by mapping arcs and isolated vertices as follows: 

(3.6) V^>cr— 1; ((z — £,j -I- £),..., (z,j)) 1-^ (i,j) and j >-^ j if j is an isolated vertex 

and secondly relabel the resulting diagram from left to right in increasing order, see Figured 
Claim 1. c: Tj^ ^(n, h) — > iJa-i<b<h-i^ki^^ — 2b,h — b) is well-defined and surjective. 
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r?T?^r\ " r?<rrT\ 



12345678 134568 123456 

Figure 4. The mapping c: Tl^i{n,h) — > Uo<6<h-iCfc("- - '2h,h - h) is obtained in 
two steps: first contraction of the stacks and secondly relabeling of the resulting diagram. 

[31 

By construction, c does not change the crossing number. Since T^ o.('^) contains only arcs of length 

[31 

> 3 we derive c{Tl. „{n)) C C1.{n — 2&, h — b). Therefore c is well-defined. It remains to show that c 

is surjective. For this purpose let 6 & C^ (n — 26, /i — b) and set a = b — {a — l){h — b) . We proceed 

constructing a fc-noncrossing structure 6 in three steps: 

Step 1. replace each label i by r^, where Vi < r^ if and only if i < s. 

Step 2. replace the leftmost arc {rp^ Vq) by the sequence of arcs 

(3.7) ((Tp - {[a - 1] + a), r, + {[a - 1] + a)), . . . , (r^, r,)) 
replace any other arc (j-p, Vg) by the sequence 

(3.8) ((Tp -[a- 1], r, + [a - 1]), . . . , (r^,, r,)) 

and each isolated vertex r^ by Tj. 

S'tep 3. Set for cc, y G Z, r;, + y < Tc + .t if and only if (b < c) or (6 = c and y < x). By construction, 

< is a linear order over 

n-2b + 2{h- b) (cr - 1) + 2a = 71 - 26 + 2{h - 6) (tr - 1) + 2(6 - (cr - l)(/i - 6)) = n 

elements, which we then label from 1 to n (left to right) in increasing order. It is straightforward 
to verify that c{6) ~ S holds. It remains to show that (5 S T^ cri^)- Suppose a contrario S contains 
an arc (i,i + 2). Since a > 2 we can then conclude that i + 1 is necessarily isolated. The arc 
(i, I + 2) is mapped by c into {j,j + 2) with isolated point j + 1, which is impossible by definition 
of C^{n', h') and Claim 1 follows. 

[31 

Labeling the h arcs of (5 £ T \. \,{n, h) from left to right and keeping track of multiplicities gives rise 

to the map 

(3.9) 

{h-b 
{o,j)l<j<h-b I X! "■?' ^ ^' Oj > O- ^ 1 



hy.T^MM^U 



.[3] 

\ It, It I r I , 

^0<b<h-l 
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given by ff, ^{5) ~ (c((5), {aj)i<j<h-b)- We can conclude that /j. ^ is well-defined and a bijection. 
We proceed computing the multiplicities of the resulting core-structures J13j : 

,r/ ^ , V^^ , 1, fb+(2-(7)(h-b)-l 

Eq. (PJU)1 and eq. ((231) imply 

(3.11) Tr,l(n,M= E ('+(',:f-^)-^)Q(n-25,/.-5), 

and the theorem follows. D 



We proceed by proving a functional identity between the bivariate generating functions of TJ^, ^(n., h) 
and C^.(n, /i). This identity is based on Theorem[3]and crucial for proving Theorem|4]in SectionlH 
Its proof is analogous to Lemma 3 in [13| . 

Lemma 2. Let fc, cr G N, fc > 2 and let u,x be indeterminants. Suppose we have 
(3.12) 

V/i>l; A,^An,h)= J2 C^^\^i^^^^^"^)^^("-26,/.-5) and A,,,(n,0) = l. 

6=cr-l ^ ^ 

T/ien we /lawe f/ie functional relation 

n>0/i<# Ti>0;i<f ^ ^ 



Proof. We set X],i>o S/i<i!. B;;(?^, h)u^x^ — S/i>o ¥'/i(^)^'' ^'^'^ compute in view of eq. (|3.12p 
(3.14) 

EEa...km-''-" = EE E B,(n-2M-fe)('+('^;;;)^(^-')-^).''.« + E-^ 

where the term Y^-^, a;* = 1-^ comes from the fact that for /i = the binomial 

5+(2-cr)(/i-6)-l 
h-b-l 
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is zero, while for any i>l the Ihs counts Ak^^in^O) = 1. We proceed by computing 

h>0b<h-ln>2h ^ ^ 



^ i:E..-.wr"^;!f;'"-')"''^^'n 



6>0 fc</i 



Setting m = h ~ h and subsequently interchanging the summation indices we arrive at 

n>Oh<-| ()>n l<m ^ ^ 



I^'^m(a;)( — 



UX^ / 1 — X 

m>0 



EEb.k^)^ 



„i'2 / 1 — a; 



X 



whence Lemma [51 D 



According to Lemma [5] and eq. p.2p we have 

(3.15) j.AnM = E r^^Tr?^"')^^^"-'''^-'^ 

(3.16) Tgi>,/.) . E r^^'rT^'^"')^^^"-^^'^-^^- 

and Lemma [5] implies the following two functional identities, which are instrumental for the proof 
of Theorem 0] in Section H) 

(3-17) eEt..(",^)"''-" = EEc.(^^,^)( "'/!tr )^" ' 



1 — a; 
n>Oh<^ 



(3.18) Et?1(«k = EEc^("''^)(t^^'^"" — 

n>0 n>0/j<# ^ 



2\cr-l ■ 

l-X 
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4. Asymptotic Enumeration 

In this Section we study the asymptotics of fc-noncrossing RNA pseudoknot structures with arc- 
length > 3 and niinimum stack length a. We are particulary interested in deriving simple formulas, 
that can be used assessing the complexity of prediction algorithms for fc-noncrossing RNA struc- 
tures. In order to state Theorem [4] below we introduce the following rational functions 

(4.1) w^{x) = 

(4.2) zo(x) = 

(4.3) uo{x) = 

(4-4) vo{x) = f \ t \2 / \ , 1 • 

X uo[x)zo{x)'' — zo{x) + 1 

Theorem 4. Let k,a €z N, k,a > 2, x be an indeterminant and pk the dominant, positive real 
singularity o/ ^^^^q /fc(2ri, 0)z" (eq. i2.5]} }. Then T[J^(n), the number of RNA structures with 
arc-length > 3 and a > 2 satisfies the following identity 

2n 



X2- 


-2 




1- 


X2 
X 




1 + 


wo(x)x^ 




Wo 




1 + 


wo{x)zo{x)'^ 


1 




zo{x) 



(A n\ Y^-r[3] f \ n / \Sr f fo n\ I ^/uo{x)zo{x) \ x zo(x) 



uoix)zo{xy - zo{x) + 1 1-x 1-zoix) 

n>0 n>0 y\/\/ \ / y v/ 



where wq{x), Z(){x),uq{x)^vq{x) are given by eq. OTll)-eg. (RTTP- Furth 



ermore 



(4.6) ^f,M-Vu{n) f^rj 

holds, where ')\.\j is the positive real dominant singularity o/ X]n>o "'"fc o-("')'^" '^'^'^ minimal real 
solution of the equation 



(. y. ^/uo{x)zo{x) 

uo(x)zo[xy - zo[x) + 1 

and tfk{n) is a polynomial overn derived from the asymptotic expression of fk{2n, 0) ^ ^k{n) ( — ) 
of eq. lOj. 



12 
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Theorem m implies the following growth rates for 3-, 4- and 5-noncrossing RNA structures with 
arc-length > 3 and stack-length > 2, 3: 



(7^'^)"' = 2.5723 



ililr' = 3.0306 



(7^'2)"' = 3.4092 



(7H)-^ 



2.0392 



(75)"^ 



2.2663 



(7^)-^ 



2.4442 



[31 [31 

Furthermore we compare in the following table the exact subexponential factors T^,^{n) (7^2)'^ 
computed via Theorem [3] with the subexponential factors ipk (n) obtained from Theorem [21 



The subexponential factor 


n 


'^f!3in)iifir 


t^U") 


Ti'U") iifir 


tZin) 


-^fkn) il'li)" 


C3(") 


50 


4.93 X 10"^ 


3.97 X lO"'' 


4.59 X 10"^ 


2.13 X 10"^° 


1.53 X 10"* 


7.81 X 10"'' 


60 


2.36 X 10"^ 


1.54 X 10"'' 


4.98 X 10"* 


3.14 X 10"" 


2.11 X 10"" 


2.93 X 10"^° 


70 


1.25 X 10"^ 


6.94 X lO"'^ 


2.86 X 10"^ 


6.22 X 10"^^ 


3.24 X 10"^° 


1.83 X 10"" 


80 


7.07 X 10"® 


3.49 X lO"'^ 


1.26 X 10"^° 


1.53 X 10"^^ 


1.71 X 10"^^ 


1.65 X 10"^^ 


90 


4.25 X 10"** 


1.91 X lO"'^ 


4.64 X 10"^^ 


4.44 X 10"^=^ 


6.00 X 10"^^ 


1.99 X 10"^^ 


100 


2.68 X 10"^ 


1.12 X lO"'^ 


1.47 X 10"^^ 


1.47 X 10"^^ 


2.98 X 10"" 


2.98 X 10"" 



Proof. In the following wc will use the notation wo,uo,zq^ eq. (|4.1[) -eq. (|4.4[) . for short with- 
out specifying the variable x. The first step consists in deriving a functional equation relating 
the bivariatc generating functions of C'^{n,h) and Chinch'). For this purpose wc use eq. p.3p 
Q(n, h) = E,<t (-l)'r?'')C,(« -3b,h- b). 
Claim 1. 



(4.8) 



5:^q(n,/^Kx" = t^EEc.K/^K(y 



n>Oh<{ 

To prove Claim 1 wc compute 



EEci(^''*) 



w'^x" 



n>n h<i 



n>0 h<\ 



EEEM)'^ .^'y.{n-ib^h-b)w^. 



EEE(-i)''('"".'^V.(--3^/^-&)-v 



EEm)"!"" '') Ec.(--3^^-^K-" 



h<^ 
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We rearrange the summation over h and arrive at 



EEM) 



n>Ofc<^ 



n - 2t> 



J2 Q(n-36,/i-6)w'' 



,.h-b 



h<\ 



W^X^ 



Setting <~Pn{w) — X]i<^ ^^.(n, j)it;-' this becomes 



E ^ — u~ E 71 ^'^"-3fc(w)a;" ^^ where to = n - 36 



6! ^^ (n-3&)! 

b>0 n>3h ^ ' 



(-wa;'^)'' .^-^ (m + fe)! 



y^ '.^"'^ ; y^ 



b>0 



6! 



m>0 



to! 



V3,„(w)x" 



'^"(^);;t E — n — ("^ + ^)! 



m>0 



m! '—' h\ 

6>0 



Laplace transformation of the series ^ 



E 

6>0 



(m + fe)! , 
6! ^ 



b>o ^^y" yields 

•^0 b>0 ' 



1 



(1 _ ,).„+! X ((l-2/)0™e-^-^M(l-y)t) 



Hence the bivariate generating function can be written as 

ipm{w)x™ ml 



^^C,*(n,/.y'x" = J2 



n>Oh<i 



m>0 



ll (1 + Wx3)™+1 



1 + wx^ ^ — ^ 1 + WX'^ 

7n>0 

T-^E E C,{n,h)w^{—^r 
1 + wic^ ^-^ ^-^ 1 + wa;'^ 



ri>0 h,<^ 
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and the Claim 1 follows. According to eq. (|3.18p . we have 



(4.9) 



EtLIh-" = EEc^("''^) T 



)i>0 



and Claim 1 provides, setting 
(4.10) 



n>Oh<i 



Wo 



(x2) 



2\cr-l 



1-X ' 



{x^Y 



the following interpretation of the rhs of eq. ()4.9I) : 



(4.11)^^Q(n,/i)f^ 



(x^Y 



1 



n>Oh< _ 

According to eq. (|3.17p and Lemma [1] we have 



1 + WqX"^ 



EECfe(">K(Y- 



n>0 h<^■ 



WqX'' 



r)" 



(4.12) ^^Tfca(n,/i)u''z" 

n>0/i<f 

(4.13) EETfc4(«>^)«''^" 



Y^Y.^kin,h)i- 



^h^n 



UZ^ 



l-z 



n>0 h<i 

and consequently 



uz^ 



-^EA-(2->0)f^^ 

- u + 1 "^-^ \ uz^ — z + 1 



n>0 



(4.14) Y: E C.(n,/.)(-^)''z" = -^^-—T E ^(2-'0) f .^^ , V" ^ T^ 
'^-^ -^-^ 1 — uz^ uz^ — z + 1 ^^-^ V uz'^ — z + 1 J 1 — z 

n>0/i<^ ri>0 ^ ^ 



n>0/i<f 

holds. Suppose now 

(4.15) 

(4.16) 



zo{x) 
uo{x) 



1 + WqX^ 
Wo 



1 + woz^ 
Then we have the formal identity "° ^ = wq and obtain substituting eq. (14. lip into eq. (14. 9p 



1 



1 + WoX'^ 



Wo 2; 



eti:l(-)-" 

n>0 ' "■" ln>0/i<f 

Next we use eq. (I4.12p and the formal identity "" ^ = wo to arrive at: 

i UnZr, 



^^C,(n,/i)«;o^- 



r 



+ 



1-a; 



EtL'U")^" 



1 



n>0 



1 + WqX'^ 



E E ^.("."KT^)'^i- 



n>0/i<-; 



l~x 
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Setting to (a;) = jj-^, ti{x) = -^ and t2{x) = -^ we derive via eq. (|4.14p 
(4.17) 

2n 



^T^S,{n)x"=toix) 



-[3] 

We compute 



,/fc(2n,0) 






-to{x)t2{x) +ti{x) . 



to{x) 



X — x^ 



x 



2CT+1 



ti{x)-toh{x) = '^ ' 



{l-x~x^+x^ + x2^+i)(a; - 1)(.t2 - 1 - x^-^+i) ' 

where /(a;) is a polynomial of degree 4cr + 2. Therefore to (2;) and ti{x) — tot2{x) are analytic in 
{zGC|0<|z|<i} and accordingly do not introduce any singularities of 

uqz^ - zq + 1 ^^ V"o2o - ^0 + 1 / 

in the domain {z£C|0<|2;|<^}. Furthermore we can conclude from eq. (|4.17p that x = is a 
removable singularity. Let us denote V{z) = X],i>o /fc(2»^,0) ( ^^;;7Tr^^j ■ 

Claim 2. All dominant singularities of X]n>o "'"L 2("')-^" ^"^^ singularities of V^(z). Furthermore the 
unique, minimal, positive, real solution of 

Uqz'^ - 20 + 1 

denoted by 7]; J^ is a dominant singularity of X)n>o 'kai^)'^^- 

Clearly, a dominant singularity of — ^^_, . y(z) is either a singularity oiV{z) or - — 13— rr. Suppose 
there exists some singularity ^ G C which is a root of ^-i^^^+i ■ ^^ construction C 7^ and ( is 
necessarily a non-finite singularity of 1^(2). If |CI 5: 7^0-1 then we arrive at the contradiction 
\V{C)\ > \V{Kk)\ since V{Q is not finite and 

7l>0 

Therefore all dominant singularities of X]n>o "'"fe cr("')'^" ^^^ singularities of V{z). According to 
Pringsheim's Theorem [2^, X]n>o "'"L cr("')-^" ^^^ ^ dominant positive real singularity which by 
construction equals 7^ J^ being the minimal positive real solution of 



(4-18) .. J . , . = Pfc 



(a;2)" 



(4-19) / (,2).-i X = Pk 

{x^Y-x^ + l 
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and the Claim 2 follows. 

[31 

Claim 2 immediately implies that the inverse oijk a equals the exponential growth-rate. According 
to [28j the power series X]n>o /'^(^'^''^)'^" ^^^ ^^^ analytic continuation in a Ap^ -domain and we 
have [z"]l/(z) ~ K(pk{n){p~^)"', where (pk{n) is given by eq. (|2.4p . We can therefore employ 
Theorem [31 which allows us to transfer the subexponential factors from the asymptotic expressions 
for //c(2n, 0) to Tk,cr{n) and eq. (|4.6p follows. This completes the proof of Theorem |4l D 
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